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ON THE FIRST SIGN CHANGE IN MERTENS’ THEOREM 


JAN BUTHE 

Abstract. The function fZ p<x ” — loglogfa;) — M is known to change sign infinitely often, 
but so far all calculated values are positive. In this paper we prove that the first sign change 
occurs well before exp(495.702833165). 


1. Introduction 

Mertens’ Theorem states that 

A M (x) := X - - loglog(x) - M = 0(log(x) -1 ) 

“ V 

p<x 

for x —>• oo, where M = 0.26149... denotes the Mertens constant [Mer74j . Rosser and 
Schoenfeld observed that A m(x) is always positive for 1 < x < 10 8 and posed the question 
whether this would always be the case |RS621 p. 72f]. This has been answered by Robin who 
showed that A m(x) changes sign infinitely often [!Rob83 ]. 

In this paper we show that the first sign change occurs before exp(495.702833165) = 
1.909875 • • • x 10 215 . More specifically, we prove 

Theorem 1.1. There exists an xo € [exp(495.702833109),exp(495.702833165)] such that 
A m{x) is negative for all x G [xq — exp(239.046541),a?o]- 

This problem is similar to bounding the Skewes number, the number in [2, oo) where the first 
sign change of A(.t) = 7r(a:) — li(x) occurs |Ske33| : this number is by now known to lie between 
10 19 [Biitl5| and exp(727.951335792) [STDarj . The functions A(x) and A m{x) are closely 
related and the Prime Number Theorem, A(x) = o(li(x)) for x —>• oo, is in fact equivalent to 
A m(x) = o(log(x) _1 ) for x —>• oo. But since A(x) and A m(x) are biased in opposite directions 
there is no correlation between the sign changes of the two functions. On the Riemann 
Hypothesis, sign changes of A m(x) rather occur at points where A(x) ~ — 2y/x/ log(x). 
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Theorem o is proven by an adaption of the Lehman method for bounding the Skewes 
number jLeh 66 j . using explicit formulas and numerical approximations to part of the zeros 
of the Riemann Zeta Function from [FKBJar] . In doing so, the kernel function in Lehman’s 
method is replaced by the Logan function |Log 88 | , which appears to be more suitable for this 
problem. This is done in such generality that it can easily be reapplied to the original Lehman 
method. 


2. Notations 


As usual C(s) denotes the Riemann zeta function and zeros of £(s) are denoted by p = (3+i'y 
with /3 ,7 € M. The Euler constant is denoted by Co = 0.57721... and the Mertens constant 
by 


OO -j 

( 2 . 1 ) " = 

to p m 

p m =2 " 

We use the symbol to define normalized summatory functions, i.e. we define 

1 \ ^ 1 

2 ^ ° n+ 2 


X' a - := \ X + \ X a - 

x<n<y x<n<y x<n<y 

Moreover, we define the Mertens prime-counting functions 


/ v ^ a * ( \ sr nM ( x / ) 

7 T M{X ) = y - and 7 t m (x) =} 

' p z —' 

p<x 1 m =1 


m 


The Fourier transform of a function / is denoted by f and defined by 

/ OO 

f(t)e~ itx dt. 

-OO 

Finally, we will use Turing’s big theta notation for explicit estimates and write f(x) = &(g(x )) 
for |/(x)| < g{x). 


3. Description of the Method 

The method we use is similar to the Lehman method for finding regions where ir(x) — li(x) 
is positive [Leh66j. We aim to calculate upper bounds for a weighted mean value 

rbJ-\-£ 

(3.1) / K(y - Lo)ye v/2 [ 7 r M (e y ) - log(y) - M] dy, 

J UJ—£ 

where K(y ) is a non-negative kernel function. By using explicit formulas this mean value 
can be expressed as a sum over the non-trivial zeros of C( s )i which can be approximated 
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K c (y) ■■= l 2 Q sinh(c) ~ 


numerically. Then, if an uj can be found for which the value in (13.11) is negative, there must 
exist an x € [exp(w — e), exp(w + e)] such that ttm( x ) ~ loglog(x) — M is negative. 

Lehman’s method uses the Gaussian function as a kernel function but we prefer to use 
dilatations of the function 

)h{c \/1 - y 2 ) \y\ < 1 , 

otherwise, 

where Io(t) := YlnL o(^/2) 2Ti /(n!) 2 denotes the 0-th modified Bessel function. The Fourier 
transform of K c is given by the Logan function (see [FKB.Tarl Proposition 4.1]) 

zv f \ _ oi \ ._ c sin(Vt 2 - c 2 ) 

Kc(x) — 4W •— • 1 r-K - - 5 

smh c V t 2 - c 2 

which satisfies an optimality property well-suited for this problem | Log88i], and which outper¬ 
forms the Gaussian function in the similar context of calculating the prime-counting function 
analytically iFKB.Jarl . 

We define 

K c ,e(y ) := - K c (y/e) and l c , £ ( x ) ■= K C)£ (x ) = £ c (sx). 

Then our main result is 

Theorem 3.1. Let 0 s < \ 10 , c ^ 3^ uj s ^ 200 ^ and let H ^ cjs be a nunxber such 

that (3 = 1/2 holds for all zeros p = (3 + of the Riemann zeta function with 0 < 7 < H. 
Furthermore, let h = 0 if the Riemann hypothesis holds and h = 1 otherwise. Then we have 

(3.2) / K Cj£ (y - a i)y e v/2 [vr M (e y ) - log(y) - M] dy 

J uj—e 

< Y. e _, 1 "4A7)fl-L 


| 7 |<c/e 


P Up* 


+ 1 + 5.4 x 10 ^ + £\ + £2 + £3, 


where 


(3.3) 

pO.Tl-^/cs / c , \ 

£ 1 <0.33e huj/2 . log(3c) logf ), 

smh c V e / 

(3.4) 

„ 3.36 + 126 e / e \«/2-i 

£zS 1000 ^ +2 ' 8 ( 2 i?) log(ff) 

and 


(3.5) 

e ^ 2 1 !tt\( ce 3 ' 12 ^ /ee\^/ 2 

£3 < QQT T l °g( H ) ■ 1 / x + ( ) 

1.99 H 1 wsmh(c) Vw/ 
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Moreover, if a E (0,1) satisfies ac/e > 10 3 in addition to the previous conditions, then 


(36) •£ 

¥<w<i 


p ujp z 


0.32 + 3.51ce /c\ cosh(c\/l — a 2 ) 

< -=-log 


car 


e/ sinh(c) 


The proof needs some preparation. 


4. The explicit formula for tt* m (x) 

The first ingredient is the explicit formula for We define the auxiliary function 


Ei(*) = [ 
Jo 


OO gZ—t 


z — t 


dt , 


which coincides with the exponential integral Ei(z) in M \ {0}, and which occurs naturally in 
explicit formulas for prime-counting functions. 


Lemma 4.1. Let x > 1. Then 

* re 

(4T) = loglog(x) + C 0 - ^2 Ei(—plogx) + / 

P Jx 
where the star indicates that the sum over zeros is calculated as 


dt 


t 2 log (t)(f 2 — 1) ’ 


lim ) Ei(—plog: 
T->oo ' 


It \<t 


Proof. The argument is similar to the original proof of the Riemann explicit formula in |vM95] . 
Let 

(4.2) 1 ° gI ’ 

Then we have 


if(x,r) = ^2 


p m <# 


P" 


7T m( x ) = J^ 4>(x,r) dr. 


From |Lan08l (39)] we get the explicit formula 

T-r _* nrp-r 


if(x, r) = 




n —2n—r 


1 — r p — r —2 n — r C 

p n= 1 


- tW 


Since Ei(— x) = log(x) + Cq + o(x) for x \ 0 [ Qlv97i p. 40], and since logC(l + e) = 
— log(e) + o(l) for e \ 0 we have 

j - - (r)dr = lim Ei(—£logx) + log£(l + e) = log log(x) + C®. 

J 1 1 r C e\o L J 
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The sum over zeros takes the form 


/ °° V |. X P ~ r V |; _ V |. 

X J—^ dr= X Ei((p-l)logx) = 22 Ei(-plogx), 


p p 

and for the sum over the trivial zeros we find 

f* OO 00 


roo 00 ^,—2 n—r 


/*oo ^ 

is 


2 n + r 


dr = 


'i 


E 

71=1 


X 


— (2n+l)r^£ _ 


/•oo 3r /-c 

L T^ ir -L 


dt 


t 2 log (t)(t 2 - 1)' 


□ 


5. The difference ir* M (x ) — ttm(x) 
By definition of the Mertens constant m we have 

= n* M (x) + M - C 0 + r M (x), 

where 


t'm(x) = Y 


/ 1 


p rn >x 
m> 2 


mp" 


The term tm{x) is responsible for the positive bias in Mertens’ Theorem and needs to be 
bounded from above. 


Lemma 5.1. Let log(x) > 200. Then 


tm(x) < 


1 + 5.3 x 10 -10 


y/x log X 

Proof. First we consider the contribution of the squares of prime numbers which yield the 
main term. Let r(t) = ip(t) — t, where if(t) := 0) in the sense of (14.21) denotes the 

normalized Chebyshov function, and assume |r(i)| < et for t > y/x and some e > 0. Then 
partial summation gives 


( 5 . 1 ) E-< 


-P 


—r(t) 1 00 

+ 


t 2 log t 


p>\/x 

For 3 < m < log(x) we use 


r_*_ r r(t )—(—-—^ 

+5 iyyt 2 log(t) J^ dt\t 2 log t) 


dt <2 


1 + 3e 
v^logx 


which gives 


E 

p>x 1 / r 

E ^ 


1 1 


< - + 


pm x 


L 


dt 1 


= - + 


„l/m—1 


l/m t m x m — 1 


p m >x 
3<m<log x 


_<!^ + ( C (2)-i )x -i/3<+54E 

mp rn x y/xiog[x) 
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For m > log x we estimate trivially: 


i 00 r 

E^sE»"”+ 2 - ”s 2 -” + 

p 1 n=3 J2 


00 rjt / 2 

= 2~ m (l + 




m — 1 


Therefore, we get 

E ^ 


1.01 v- m 2.02 x 2 _log ( a: ) 1(T 16 

< - \ 2 < --——- < - 

mp m log(x) 

p m >x y 6V ' m>log x 

m>log x 


log(x) v / xlog(x) 

By [ Biitar , Table 1] ()5.1D holds with e = 1.752 x 10 -10 and so the assertion follows. 


□ 


6. Evaluating the sum over zeros 

The next problem is to approximate the following integral of the sum over zeros 

K c ,e{y - u)ye y/2 ^2 Ei (-py)dy. 


j: 


Here, integral and sum may be interchanged, since the sum converges locally in L 1 . Therefore, 
we may treat each summand individually. 

6.1. Asymptotic Expansion of the Summands. Since the Logan kernel should also be of 
interest for the question on finding regions where n(x) — li(x) is positive, the following Lemma 
is presented in a more general version, which also covers the classical case. 

Lemma 6.1. Let 0 < s < ui, and let K € L 1 ([— e,e\) satisfy ||it || L i = 1. Let a € [0,1], let 
p = /3 + i'y, where 0 < (3 < 1 and 7 € M \ {0}, and let 

fUl+E 


Then we have 

( 6 . 1 ) 


r^-\-£ 

$u, P ,a= / K(y-uj)ye ( 2 ~ a)y Ei((a-p)y)dy. 
J LJ—£ 


FuJ\o) q ( k\e 2 


$u, P ,a 1)! ( p - a )j 

j=i 


(a; - e) k |y| 


k U| fe +1 


where F^p^if)) = — and for j > 2 and any m > 0 , 


(6.2) fW>(0) = (-1)A<U<’)»^ 


71=0 


n 


+ j-2\ (-»)"£<")(f-E) 


n 


u; n +j 


-1 


+ 0 


e? 2+ 2 e ^2 (ee/u) 


771+1 




-1 


1 — ee/uj 
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Proof. By definition of Ei we have 

ruj+e 


oo e (a-p-r)y 

a — p — r 


dr dy 


(6.3) 

Now let 


ruj-\-e r 

*, P ,a= K{y-u)ye^y 

J LO — £ J 0 

r oo i rw+£ 

= / - / K(y — dydr. 

Jo a- p-r J u _ £ 

ru+e 

Fpp(r) := (-1 y / y l - j K{y - U )e^~^ dy, 

J U)—£ 


which is well defined since u > e, and satisfies pF^p = Fp p ^. Then partial summation 


gives 


$ 


u,p,a 


Jo r + p-a 
Here, the trivial bound 




3 = 1 


(p ~ a) j 


f 


dr. 


(r + p- a) k+1 


Fp p k) (r) | < J 


c (b-0~r)(y+u) 


-e (^ + v) k 1 


dy < 


£f_ e (h-P)u e r{£-u) 


(uj — e) k ~ l 


yields 


poo 

F^p k \r) 

^ eze^2 P) u 

L 

(r + p- a) k+1 

' ^ (ui — e) k |7| fc+1 


which confirmes (16.11) . It remains to evaluate F^J\o). For j = 1 we find 
Fpp (0) = J K(y)e~ i ^-^y dy = -e^~ p ^K^ - 

For larger values of j we use the Taylor series expansion 


1 


(w + y) 1 


E 

n=0 


u + n - 1\ (-y) n 
n ) u u+n 


and 

(6.4) 

which gives 


f K(y)y n e~ i ^-^y dy = i n K (">(£ - 1), 

J + n-2^ (-*)"£(">(?-£) 

-l 


^ ) ( 0 ) = (-i) ie(3 “ p)w E 

n=0 

From (16.41) we get 

IfWf^ - -) 
\i 2iJ 


n 


u n +J 


< ei£ r ‘ 
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and the inequality (^) < ) b , which follows from Stirling’s lower bound for 6!, implies 


i+n - 2 W(i + ^'i < e -+'- 2 . 

n V n 


Thus, we have 

OO 

£ 

n=m -\-1 


j + n — 2 
n 


A <n) (f - s) 


which confirms the bound in (16.21) . 


< e - 7 2+ 2 ^ /ee\ n e J 2+ 2 (ee/cu ) m+1 

— V u) a;-? -1 1 — ee/uj ’ 

n=m+l 


□ 


6.2. Bounds for the Kernel Function. We need some bounds to estimate the tails of the 
sum over zeros. These are provided by the following two Lemmas from [Blit] and [Biitarj : 


Lemma 6.2 ( [Biitarl Lemma 2]). Let 0 < e < 10 3 and c > 3. Then we have 
(6.5) 


£ . 2i)l < 0.32fLLl log(3c) log( 

It! smh(c) v 


Lemma 6.3 ( [But l Lemma 4.5]). Let 0 < e < 10 3 and c > 3, and let a E (0,1) satisfy 
ac/e > 10 3 . Then we have 


( 6 . 6 ) 


£ 


\Ic,z{t)\ < 1 + lice 1 fc ^ cosh(c\/l — a 2 ) 




ItI 


7rear 


log 


sinh(c) 


We also need bounds for the derivatives 1™) (^ — ?p) occurring in (16.21) . for calculations not 
assuming the Riemann hypothesis. 

Lemma 6.4. Let 0 < e < 5 < c/100, and let z € C satisfy |3?(^)| > c/e and |9(z)| < 4. Then 




< n! 


ce 


1.56 y/Sc 


sinh(c) V 5 

Proof. The bound follows from the Cauchy formula 


/2e\ n 


4*2 (*) = — / 

C > £ V 2m 


r(0 


i*-*i=£ (* - ^) n+1 


dt, 


if we show that 
(6.7) 


IMOI < 


ce 


1.56\/<5c 


sinh(c) 
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in the range of integration. By basic properties of t Cj£ it suffices to prove this bound for e = 1 
under the conditions 9i(£) > c — 5, 0 < SS(£) < 5, and we may also assume 5 < c/100. Since 
we have 


9?(-\/£ 2 — c 2 ) < 9(y/(c — b + ib) 2 — c 2 ) 

/—7--t-t- /7T 1 (b C—5 2 

< v/2 1 + i\bc sm — H— arctan —-- 

~ V 4 2 V <5c 

< 2 3 / 4 sin(1.181)\/^c < 1.56>/£c 


under these conditions, the desired bound follows from 


sin(z) 


< e 


|9«l 


□ 


7. Proof of Theorem 13.11 
By Lemma 14.11 and Lemma 1-5. II we have 

KM(e y ) ~ log (y) — M = Tr* M (e y ) - log (y) - C 0 + r M (e y ) 

, 1 + 5.4 x 10” 10 „ /2 

<-V Ei {-py) + - e~ y ' 2 

p y 

for y > 200, where we estimated the integral in (14.11) trivially by e~ 3y . Therefore 

ruj-\-£ 

/ K c>£ (y - u)ye y/2 [ir M (e y ) - log(y) - M] dy < - ^ <L w , P ,o + 1 + 5.4 x 10~ 10 , 

J UJ—£ 


with $ WjP) o as defined in Lemma 1671 with K = K Cj£ and K = i C}£ . We subdivide the sum over 
zeros into two parts. For 0 < 7 < H we choose k = 2 and m = 0 in Lemma 16.11 which gives 


( 7 . 1 ) E 

\l\<H M<c/e f<H < H 


4,6 ( 7) 


7 


1 + — 

CUJ 


+ 772 E 


2.72e 2.01 

—o—b 


l7l<« 


br 
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where we have used e < 10 . For 7 > H we have 

k 

(7.2) Y, l^, P ,ol <e hw/2 Y ' 


l7l>« 

k 

+e hu/2 y 

h\>Hj=2 


|7l >H 


(J- 1 )! 


V z 2z > 


7 




1=1 


uj 3 


-1 


fcS ItP 


E 

1 n =1 


fn + j- 2\ 

M)ip _ 1 1 
2 i> 

l n I 

u n + 3 ~ l 


ra +1 


+ 


+ e 


e j 2+e/2^ e£ 'j 

u 3+m - 1 {u - ee) 

k\e £ / 2 


hcu/2 


E 

l7l>« 


(w-e) fc |7 | fc+1 


for arbitrary k > 2 and m > 1, where h = 0 if the Riemann hypothesis holds and h = 1 
otherwise. So the inequality in (13.2 jl holds with 


(7.3) 

(7.4) 
and 


*= E 

§<l7l<ff 


4 +(7) 


i +—) +e huj/2 y 

cwJ , 4 ^ 

l7l>-ff 


4,e(i 2i) 


E 


U ~ !) ! Tjl- 


3 = 1 


+4 


-1 


H 3 , 


£2 = 


1 ^/2.72e 2.01 


OJ 


2 Z-/V <1/2 


i7r 


+ e' 


hcj/2 


E 


k\e £ / 2 


(w-e) fc |7| 


fcUl fc+r 


(7.5) £ s = e“/ 2 ^ 53 

|7l>^J = 2 


(J-1)! 

fcS I7P 


E 

1 n=1 


n + j — 2 
n 


p(n) ( p Js 
<7+ l j 2t / 


of —2 + e , 


W n +J 


-1 


+ 


/v: 


m+l 


u;J+ m - 1 (w - ee) / ' 


We proceed by bounding £*.. To this end we choose k = m = |u>/2j. In (17.31) we take 
H = £ which gives 


(7.6) 


4,6 (7) 


7 


fe-i 


£\ < e ha; / 2 

I <ItI 

where the inner sum is bounded by 

E4)E(i- 

1=0 

since c > 3. Using this and (|6.5I) in (17.61) gives (13.31) . 


n I 


E^(), 

1=0 


1 N " 1 


6000 


< 1 . 0002 , 


In (17.41) we use the bounds ^ 7 2 < 0.0463 and ^ I 7 I 3 < 0.00167 from |Ros411 Lemma 


17], the bound 


E br* <^logcr) 

l7l >T 


(7.7) 
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for T > 2ire and k > 2 from [Leh 66 l Lemma 2], and the inequality (u — e) k > e £ u k , which 
follows from k < w/ 2 , and get 


£2 < 

In (1731) 

gives 


0.00336 + 0.126e 




+ e‘ 


w.jfIE log(ff) <?4™ +2 , 




1000 a ; 2 


(£) 


ui/2-1 


log(H). 


we use (17.71) again and the bound from Lemma 16.41 where we choose <5 = 4e. which 


(7.8) 


J=2 




sinh(c) “—' cu 

v ' n= 1 




+ 


)-2 


1 . 002 e ^~ 1 (j - 1 )! /e£\ m + 1 
e c4 _1 V w / 


Since n+j — 2 < u we have (n+j — 2)\/u n+ ^~ 2 < 1 /w, so the inner sum is bounded by l/(2cu). 
In the second summand, we use the bound (j — 1)!/co - 7 — 1 < 2 1- - 7 . Since Yl'jLi H ^ — 1-001/ H, 
2H/e)-i < 1.001e/(2iL), and m + 1 > u/2, we obtain the bound in (13.51) . 

Finally, the estimate in (13.61) follows from (16.61) since 


E 

f<l7l<| 


4 , 47 ) ^ _ J_\ 

p V up) 


< 



1 


200 x 1000 


E 

f <l7l<f 


4 , 47 ) 

7 


□ 


8. Numerical Results 

To locate potential regions where the left hand side of (13.21) should be small, the function 

e +y 

■ — i'y 

has been evaluated for T = 10 6 at all points in 10 -7 Zn [1,2500]. Since 4, 5 ( 7 ) = l+0((e7) 2 /c) 
for 7 —>■ 0 this gives a reasonably good approximation to the first part of the sum in (13.2p . 
and the objective is thus to find regions where 07 , (y) smaller than — 1 . 

The evaluation has been done using the method for fast multiple evaluation of trigonometric 
sums from IFKB.Iar I. A more detailed search with T = 10 8 around 495.7028078, the first point 
where cr W 6 ( y ) turned out to be promisingly small, revealed a short region of length ~ 2.8 x 10 -8 
about 495.702833137 where cr w s(y) is smaller than —1. 

Theorem 11.11 now follows by an application of Theorem 13.11 with u = 495.702833137, c = 
280, e = 2.8 x 10~ 8 , H = 10 11 (which has been reported in [FKB.Iarj l and a = 0.4. 


°t{v) = 
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Table 1. Values of y € [1,2500] for which <r 10 6 (y) < —0.95. 


_ y_ 

495.7028078 

1423.957207 

1623.9204309 

1859.1291846 

2107.5263606 

2285.3917834 

2430.3039554 

2447.6661764 


Q~iq6 (y) 

-0.9972... 
-0.9740... 
-0.9807... 
-1.0511... 
-1.0214... 
-1.0454... 
-1.0172... 
-1.0028... 


The sum over zeros was calculated using approximations to the zeros with imaginary part 
up to 4 x 10 9 which were given within an absolute accuracy of 2 -64 . The sum was evaluated 
using multiple precision arithmetic, which gave the bound 


(8.1) 


|7|<4x10 9 ^ 


— 7 < -1.00015419. 

cop z ) 


The sum in (13.61) is then bounded by 1.2 x 10 11 and we have 

Ex + £2 + £3 < 1.2 x 10“ 12 + 1.37 x 10~ 8 + 1.6 x 10~ 24 < 1.38 x 10" 8 . 


Thus, the left hand side of (13.21) is bounded by 

-1.00015419 + 1.2 x 10" 11 + 1 + 5.4 x 10" 10 + 1.38 x 10" 8 < -0.000154. 


Consequently, there exists an x E [exp(tc — e), exp (w + e)] such that 7 tm{x) — log log(x) — M < 
—0.000154/(y / xlogx). Obviously, we have 


km(x - y) - log log(x — y) — M < ir M (x) - log log(x) 



dt 

t log t 


0.000154 y 

<- 1 - - - 5 

v / xlog(x) (x - y) log(a; - y) ’ 

which is negative for y < 0.00015\/T. Since O.OOOlS^/x > exp(239.046541) the assertion of 
Theorem o follows. □ 
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